Mott— Hubbard metal— insulator transition at non-integer filling 
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Correlated electrons in a binary alloy A^Bi-^ are investigated within the Hubbard model and 
dynamical mean-field theory (DMFT). The random energies have a bimodal probability distribu- 
tion and an energy separation A. We solve the DMFT equations by the numerical renormalization 
group method at zero temperature, and calculate the spectral density as a function of disorder 
strength A and interaction U at different fillings. For filling factors v = x oi 1 -\- x the lower or 
upper alloy subband is half filled and the system becomes a Mott insulator at strong interactions, 
with a correlation gap at the Fermi level. At the metal-insulator transition hysteresis is observed. 
We also analyze the effective theory in the A — > oo limit and find good agreement between analytical 
and numerical results for the critical interaction Uc at which the metal-insulator transition occurs. 

PACS numbers; 71.10.Fd, 71.27.+a, 71.30.+h 



I. INTRODUCTION 

A Mott-Hubbard metal-insulator transition (MIT) oc- 
curs when the local interaction U between the electrons 
on a lattice reaches a critical value Uc<^ Since in the 
Mott insulator the particles are essentially localized on 
lattice sites, the transition can only occur if the number 
of particles is commensurate with the number of sites, 
i.e., for integer filling factor v = N/Nl, where N and 
Nl are the number of electrons and lattice sites, respec- 
tively. If the translational symmetry of the lattice is bro- 
ken either spontaneously (as in an antiferromagnet) or by 
application of some field (as in a superlattice) the unit 
cell is enlarged and the Mott insulator can, in principle, 
occur even at a rational filling factor v — p/q, where 
p,q G NiSi^i^ In fact, in our recent investigation of dis- 
ordered electronic systems^ we discovered that in binary 
alloy systems A^Bi-x, composed of two different atoms 
A and B with concentrations x and 1 — x, respectively, 
a Mott-Hubbard MIT can even occur for arbitrary non- 
integer filling factors. This transition takes place if — x 
or 1 -|- x, where < a; < 1. This observation extends the 
traditional view of the Mott-Hubbard MIT and the no- 
tion of Mott insulators to a wider class of systems with, 
basically, arbitrary fillings. As we will argue below, in 
such a system magnetic long-range order should be sup- 
pressed due to the absence of particle-hole symmetry 
and, therefore, the ground state is likely to be param- 
agnetic. Therefore, the experimental realization of such 
a Mott insulator would provide an excellent playground 
for the study of Mott insulators without long-range or- 
der. 

In the present paper we study correlated electrons 
on a lattice using the dynamical mean-field theory 
(DMFT)'^''^ applied to the disordered Hubbard model 
at zero temperature^ The numerical renormalization 
group (NRG) method is used to solve the selfconsis- 



tent DMFT equationsAifliiiiiS We present results for the 
single-electron spectral density and the self-energy to 
show how a correlated metal away from half-filling may 
turn into a Mott insulator by increasing the disorder. We 
find hysteresis in this transition, like the one observed for 
the MIT in the pure Hubbard model at half-filling 
From the numerical data we construct a phase diagram 
at v — X and propose a new classification scheme for 
correlated insulators with binary-alloy disorder. We also 
develop an effective analytic theory which is asymptoti- 
cally exact in the limit of alloy band splitting, and which 
shows that the Mott-Hubbard MIT occurs at the critical 
interaction Uc/W « 3y/x/2, where W is the band-width 
of the noninteracting system. 

In Ref. we analyzed the problem of the Mott MIT 
in a binary-alloy host at relatively high temperatures us- 
ing the quantum Monte-Carlo (QMC) method to solve 
the DMFT equations. This only allowed us to detect a 
crossover from a correlated metal to a Mott insulator. 
Moreover, since in Ref. |^ we were primarily interested 
in ferromagnetism in binary-alloy systems, we used the 
density of states (DOS) corresponding to the fee lattice 
in infinite dimension. Such an unbounded DOS sup- 
ports ferromagnetism in a one-band Hubbard modeliSiii 
but does not lead to a real gap for the Mott insulator. 
The NRG method applied here is accurate enough to de- 
tect a sharp MIT in a disordered and correlated electron 
system at zero temperature. Explicit calculations were 
performed with a semi-elliptic Bethe DOS, having finite 
support, which leads to the opening of a genuine gap at 
the MIT. Additionally, the NRG allows us to determine 
the self-energy on the real frequency axis and to discuss 
in detail how the Mott gap opens. 

The paper is organized as follows: in section II we 
introduce the model and discuss the properties of inter- 
acting electrons in a binary alloy host. In section III the 
DMFT method is introduced to solve the Hubbard model 
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with local disorder. In section IV, we present the numer- 
ical results and provide evidence for a MIT away from 
half-filling. In section V an asymptotic theory is devel- 
oped, and in section VI our conclusions are presented. 
In the Appendix, we present a novel method to extract 
the local self-energy within the DMFT in disordered sys- 
tems. 



II. ANDERSON-HUBBARD HAMILTONIAN 
WITH BINARY ALLOY DISORDER 

As a minimal model describing correlated lattice elec- 
trons in the presence of disorder we consider the single- 
orbital Anderson-Hubbard (AH) Hamiltonian 

Hah = "i<T«j'7 + X! + ^ X!"*T"a' (1) 

{ij)cr ia i 

where t > is the hopping integral for the electrons be- 
tween nearest neighbor sites, U is the on-site interaction 
energy between electrons with opposite spins a = ±1/2, 
nia = ^hc local electron number operator, and 

Si is the local ionic energy which is a random variable. 
In the following we assume a bimodal probability distri- 
bution for the random variable e^, i.e., 

V(e,)^xS{e. + ^) + {l-x)S{e,-^), (2) 

which corresponds to a binary-alloy system composed of 
two different atoms A and B. The atoms are distributed 
randomly on the lattice and have ionic energies ea.b, 
with — Eyi = A. The parameter A is a measure of the 
disorder strength. The concentration of A (B) atoms is 
given by X = Na/Nl (1 - x = Nb/Nl), where Na (Nb) 
is the number of the corresponding atoms. 

From the localization theorem (the Hadamard- 
Gerschgorin theorem in matrix algebra) it is known that 
if the Hamiltonian Hah, with a bimodal distribution for 
€{, is bounded, then there is a gap in the single-particle 
spectrum for sufRciently large A ^ max(|t|, U)ii^^ 
Hence at A = Ac the DOS splits into two parts cor- 
responding to the lower and the upper alloy subbands 
with centers of mass at the ionic energies and es, re- 
spectively. The width of the alloy gap is of the order of 
A. The lower and upper alloy subband contains 2xNl 
and 2(1 — x)Nl states, respectively. If the Hamiltonian 
is not bounded, as for example in the case of a tight 
binding Hamiltonian on a hypercubic lattice in infinite 
dimensions, the alloy gap is reduced to a pseudo-gap, 
i.e., the spectral function vanishes only at a single point. 
These statements hold for all space dimensions. However, 
the alloy gap can be destroyed by clusters of one type of 
atoms which are surrounded by atoms of the other type. 
They create an exponentially vanishing DOS in the gap 
(Lifshitz tails) near the edges of the alloy subbandsii^ 

Binary alloy disorder in a noninteracting electron sys- 
tem can create two kinds of localized states: (i) states 



which are localized due to coherent backscattering pro- 
cesses (Anderson localization) jSSiSi and (ii) states in the 
middle of the alloy subbands, which are localized due to 
a particular superposition of the electronic wave func- 
tions caused by particular arrangements of the alloy 
atomsiiSiSSi^ While the localized states of type (i) are 
generic and gradually appear in the alloy subbands start- 
ing from the band edges, the localized states of type (ii) 
can be removed either by a small perturbation of the ionic 
energies or by an interaction between the particle, and, 
therefore, are beyond the scope of the present paper. 

The most spectacular effect of strong correlations be- 
tween the electrons in a pure (A = 0) system is the 
Mott-Hubbard metal-insulator transition. It occurs for 
i^ = l,2,...,25 — 1, where g is an orbital degeneracy, at 
an interaction U = Uc (note that v = 2g corresponds to 
a band insulator) ii2i2S*2SiSL2S For interactions U Uc 
the electrons gain kinetic energy due to the delocalized 
nature of the wave function. When U increases, the 
electrons keep at a distance as much as possible, and 
'At U > Uc the many-body wave function is essentially 
localized at each lattice site. The system is then a Mott 
insulator. Within the DMFT for Hubbard-like models 
with single or degenerate orbitals a Mott-Hubbard MIT 
was shown to occur at integer filling factors i/iSL2& At 
low but finite temperatures this MIT is discontinuous, 
while at zero temperature it is continuous, i.e., when ap- 
proaching the critical interaction from the metallic side 
the quasiparticle peak continuously narrows until it com- 
pletely disappears at the transition point ii^ii In the 
insulating phase the DOS is zero at the Fermi level and 
the whole spectral weight is shifted into the two Hubbard 
subbands, which are remnants of the atomic levels with 
single and double (multiple) occupancy. Away from inte- 
ger filling the system is always metallic and the spectral 
function has visible Hubbard subbands at any finite and 
large U. 

New possibilities appear in systems with correlated 
electrons and binary alloy disorder. The Mott-Hubbard 
metal insulator transition can occur at any filling v ~ x 
or 1 + X, corresponding to a half-filled lower or to a 
half-filled upper alloy subband, respectively, as shown 
schematically for v — x in Fig. The Mott insula- 
tor can then be approached cither by increasing U when 
A > Ac (alloy band splitting limit), or by increasing A 
when U > Uc (Hubbard band splitting limit). The na- 
ture of the Mott insulator in the binary alloy system can 
be understood physically as follows. Due to the high en- 
ergy cost of the order of U the randomly distributed ions 
with lower (higher) local energies are singly occupied 
at v = X (v ~ I + x), i.e., the double occupancy is sup- 
pressed. In the Mott insulator with v — x the ions with 
higher local energies are empty and do not contribute to 
the low-energy processes in the system. Likewise, in the 
Mott insulator with v = 1 + x the ions with lower local 
energies are double occupied implying that the lower al- 
loy subband is blocked and does not play any role. We 
note that for finite U virtual processes leading to double 
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FIG. 1: Schematic plot representing the Mott-Hubbard 
metal-insulator transition in a correlated electron system 
with the binary alloy disorder. The shapes of spectral func- 
tions A{u}) are shown for different interactions U and disorder 
strengths A. Increasing A at f/ = leads to splitting of the 
spectral function into the lower (LAB) and the upper (UAB) 
alloy subbands, which contain 2xNl and 2(1 — x)Nl states 
respectively. Increasing 17 at A = leads to the occurrence 
of lower (LHB) and upper (UHB) Hubbard subbands. The 
Fermi energy for filling v = x is indicated by /i. At = x 
(or V = 1 -\- X, not shown in the plot) the LAB (UAB) is 
half-filled. In this case an increase of U and A leads to the 
opening of a correlation gap at the Fermi level and the system 
becomes a Mott insulator. 



occupation either in the lower [v = x) or in the up- 
per alloy subband {v = 1 -\- x) are still possible, lead- 
ing to the antiferromagnetic superexchange interaction. 
However, since the positions of the corresponding atoms 
are random, particle-hole symmetry is absent such that 
long-range antiferromagnetic order is expected to be sup- 
pressed. Of course, antiferromagnetism cannot be ruled 
out completely on these reasoning and its appearance in 
the model requires further studies. 

For U > J7c(A) in the Mott insulating state with bi- 
nary alloy disorder one may use the lowest excitation 
energies to distinguish two different types of insulators. 
Namely, for U < A an excitation must overcome the 
energy gap between the lower and the upper Hubbard 
subbands, as indicated in Fig. We call this insulat- 
ing state an alloy Mott insulator. On the other hand, for 
A < [/ an excitation must overcome the energy gap be- 
tween the lower Hubbard subband and the upper alloy- 
subband, as shown in Fig. ||2Jl. We call this insulating 
state an alloy charge transfer insulator. 



III. DYNAMICAL MEAN-FIELD THEORY FOR 
THE DISORDERED HUBBARD MODEL 

The Mott-Hubbard MIT is driven by the interaction 
between the particles. Since this transition occurs when 
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FIG. 2: Two possible insulating states in the correlated elec- 
tron system with binary-alloy disorder. When [/ < A the 
insulating state is an alloy Mott insulator with an excitation 
gap in the spectrum of the order of U . When 17 > A the 
insulating state is an alloy charge transfer insulator with an 
excitation gap of the order of A. 



the interaction energy is comparable with the single- 
particle energy of the electrons, there is no natural small 
parameter (e.g. t or U) in the theory. The problem is 
generically nonperturbative. Moreover, when the transi- 
tion appears between a paramagnetic metal and a para- 
magnetic insulator, there is no obvious order parameter 
characterizing the insulating phase. In the following the 
insulating state at T = is defined by the vanishing of 
the one-particle spectral function at the Fermi level. 

The model defined by the Hamiltonian |^ is not ex- 
actly solvable for any finite number of space dimensions. 
However, with a proper rescaling of the hopping integral 
it becomes numerically solvable in infinite dimensions, 
i.e., within DMFTiSSiSii For finite-dimensional systems 
DMFT is a self-consistent approximation scheme which 
takes into account local quantum fluctuations but ne- 
glects spatial correlations. Since DMFT is a nonpertur- 
bative method, it is ideally suited to study the Mott- 
Hubbard MIT. 

To derive the DMFT equations for the problem at hand 
we select a single lattice site, i, and integrate out all 
the electronic degrees of freedom corresponding to other 
sites.- This leads to an effective single-impurity Ander- 
son Hamiltonian 



H 



SIAM 



= ^(ti - lJ^)a\„ai„ + Unipi,i (3) 



ktr 



kCT 



where fj, is the chemical potential, Vk and ek are the 
hybridization matrix element and the dispersion rela- 
tion for the auxiliary bath fermions Cko- , respectively. In 
the present paper the Hamiltonian is solved at zero 
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temperature using the numerical renormalization group 
methodiSiiSiiiii^ For each ionic energy we obtain the 
local Green function G(aj,ei). The physical Green func- 
tion ^ is obtained by algebraic averaging of G{uj,ei) 
over different realizations of the disorder^ i.e., 



G{uj) = I de,r{e,)G{iu,e,). 

From the k-integrated Dyson equation 
G'-^(w) cj - 77(cj) - 



(4) 



(5) 



we determine the local self-energy The function 

r]{uj), called hybridization function, is defined as 



(6) 



The DMFT equations are closed by a Hilbert transform, 
relating the local Green function for a given crystallo- 
graphic lattice to the self-energy, i.e.. 



G{uj) = / de 



Noie) 



(7) 



where iVo(e) is the non-interacting DOS. Eqs. (0HZIl con- 
stitute a closed set of DMFT equations for the disordered 
Hubbard model. 

We solve Eqs. for a Bethe lattice with infinite 

connectivity. In this case the DOS is given by 



(8) 



with the bandwidth W — 2D. In the following we set 
D = 1. With the DOS I© the Hilbert transform (0 can 
be calculated analytically leading to a simple algebraic 
relation between the local Green function G(lo) and the 
hybridization function 77(w), i.e.. 



(9) 



In the numerical calculations we adjust the value of 
the chemical potential \x so that the number of particles 
in the system is fixed. Hence the independent variables 
are z^, ?7, x, and A. 

Since DMFT neglects short-range spatial correlations, 
and hence does not include effects due to backscattering 
of electrons on the randomly distributed ions, it cannot 
describe effects of Anderson localization. On the other 
hand, binary alloy disorder is a particularly strong type 
of disorder since it even leads to band splitting - and 
thereby to insulating behavior - in any spatial dimen- 
sion. This dominant feature, and all other effects caused 
by the simultaneous presence of interactions and disor- 
der, is well captured by DMFT. In particular, the DMFT 
equations reduce to the equations of the coherent poten- 
tial approximation (CPA) for interaction U = 0»2& The 



CPA method is known to be very successful in explain- 
ing single-particle properties of disordered systems, both 
in the case of models and realistic systemsii2i^ In par- 
ticular, it reproduces the alloy band splitting in binary 
alloy systems. Therefore we use the DMFT to describe 
the Mott-Hubbard MIT in the presence of binary alloy 
disorder, and consider additional effects due to Anderson 
localization as secondary. 



IV. NUMERICAL RESULTS FOR THE 
DISORDERED HUBBARD MODEL 

In the following we present our numerical results for 
filling factors v ^ x and v ~ x with equal concentra- 
tion of A and B atoms, i.e., x = 1/2. In particular, the 
ground-state phase diagram at = a; is presented and 
the MIT is discussed in details. 



A. Interacting electrons in tlie alloy band splitting 
limit for v ^ x 

The influence of the disorder-induced alloy band split- 
ting on the spectral function A(lo) is shown in Fig. |(2Jl for 
U ~ ^ and filling factor v ~ 0.3. For vanishing disorder 
(A ~ 0) the spectral function is composed of the lower 
Hubbard subband (which at this low density is merged 
with the quasiparticle peak) and the upper Hubbard sub- 
band around uj — A. Upon increasing A the upper alloy 
subband splits off and moves to larger ui. At the same 
time the lower alloy subband appears with a smaller num- 
ber of states. Its shape and the position with respect to 
the Fermi energy does not change for A > 2. These re- 
sults and the shape invariance of the lower alloy subband 
in the presence of interactions between the electrons sug- 
gest that in the A — > oo limit the binary alloy disordered 
Hubbard model can be effectively described by the Hub- 
bard model in a reduced Hilbert space which contains 
only the states from the lower alloy subband. 



B. Phase diagram and Mott-Hubbard transition at 
filling = X 

By solving numerically the DMFT equations we ex- 
tracted the zero temperature phase diagram at = a; = 
0.5 which is presented in Fig.^ The curve with filled dots 
represents the critical interaction Uc2 = C^c2(A) separat- 
ing the paramagnetic metal (PM) and the paramagnetic 
insulator (PI). This boundary line was determined by 
solving iteratively the DMFT equations using a metallic 
hybridization function as an initial input. It means that 
in solving the system of Eqs. iteratively we began 

with the hybridization function ■q'^^^cj) « G^^^uj) that 
has a nonvanishing imaginary part at a; = 0. The other 
curve (open dots) represents the boundary Ud = Uci{A) 
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FIG. 3: Spectral function of the Hubbard model ai U — 4 
and v = 0.3 for different binary alloy disorder strengths A. 
The upper alloy subband splits off for A > 2. The shape of 
the lower alloy subband is not significantly changed for large 
A. 



between the metallic and insulating phases, as deter- 
mined by solving the DMFT equations with an insulating 
hybridization function as an initial input. In this case 
the initial hybridization function had vanishing imagi- 
nary part at w = 0. The boundary points (lines) corre- 
spond to the values of the {U, A) parameters where the 
converged spectral functions obtained from Eqs. m?)) 
starts to have zero weight at w = 0. In the inset to Fig. 01 
the behavior of the spectral function at the Fermi level is 
shown when the metal-insulator boundary is approached 
from the metallic (solid lines) and from the insulating 
side (dashed lines), respectively. We observe hysteresis, 
which indicates that at low but finite temperatures the 
MIT transition in the Hubbard model with binary alloy 
disorder is discontinuous. By contrast, at zero temper- 
ature the transition should be continuous and occurring 
at the Uc — ?7c2(A) lineii^ We calculated numerically 
the average density of double occupied sites d = (n^ni) 
in the coexistence regime and found that d is larger for 
a metallic solution. It implies that the metallic ground 
state is energetically more stable.^ Of course, very close 
to Uc2iA) we cannot make an absolutely precise state- 
ment because of the finite numerical accuracy. 

From the inset to Fig. ^ we also conclude that in the 
metallic phase A{0), the spectral function at the Fermi 
level, decreases with disorder but remains independent 
of U. This behavior corresponds to the "pinning" of 
the spectral function (Friedel sum rule) in the pure case, 
where A{0) does not depend on the interaction U be- 
tween the electronsi^ Similar behavior is encountered in 
the disordered Hubbard model studied here. However, 
now A(0) is reduced by the A dependent imaginary part 
of the self-energy. 

In the upper panels of Figs. I5I7I we present the spec- 
tral functions for selected parameters U and A along a 
path in the ([/, A) phase diagram (Fig. ^ indicated by 



FIG. 4: Ground state phase diagram of the Hubbard model 
with binary-alloy disorder at filling z/ = x = 0.5. The filled 
(open) dots represent the boundary between paramagnetic 
metallic (PM) and paramagnetic insulating (PI) phases as de- 
termined by DMFT with the initial input given by the metal- 
lic (insulating) hybridization function. The horizontal dotted 
line represents Uc obtained analytically from an asymptotic 
theory in the limit A — > oo (see section V). Stars show the 
points at which the explicit spectral functions are presented 
in Figs. |21-|7| Inset; hysteresis in the spectral functions at 
the Fermi level obtained from DMFT with an initial metal- 
lic (insulating) host represented by filled (open) symbols and 
solid (dashed) lines. 



crosses: [/ = 3 and A = ^ 1 ^ 1.5 -> 2 (Fig. E|, 
A = 2andC/ = 3^2 (Fig. EJ, and finally, C/ = 2 and 
A = 2 ^ 1 ^ (Fig.d). These spectral functions illus- 
trate the evolution of the system within, or between, a 
metallic and an insulating phase when disorder and in- 
teraction are changed. In particular, in Fig. we see 
for U = 3 how increasing A from the value (where 
the spectral function is qualitatively similar to that of 
the v = 0.3 case, c.f. Fig. to 1 and 1.5 leads to the 
splitting of the alloy subbands and the emergence of the 
quasiparticle peak at w = 0, a feature of strong correla- 
tions between the electrons. At A = 2 the quasiparticle 
peak is absent and the spectral function possesses a Mott 
gap at the Fermi level, a feature of an insulator. Keeping 
A = 2 in Fig. El and lowering U from 3 to 2 leads to a 
shrinking of the Mott gap and reappearance of a quasi- 
particle peak, characterizing a correlated metallic phase. 
Finally, upon lowering A from 2 to at constant U — 2 
the alloy subbands approach each other and the quasi- 
particle peak merges with the lower Hubbard subband as 
presented in Fig. |7| 

In addition to the spectral functions, the imaginary 
and real parts of the self-energy, calculated by the 
method presented in the Appendix, are shown in pan- 
els (b) and (c) in Figs. EHZI In the metallic phase with 
A > 0, two important features of the self-energy should 
be noted: (i) the imaginary part of the self-energy at 
the Fermi level is finite, i.e., ImE(a' = 0) < 0, and (ii) 
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the real part of the self-energy at the Fermi level has a 
negative slope, i.e., dRe'E{uj = 0)/duj < 0. While the 
former feature (caused by disorder even in the presence 
of the local interaction) can be observed within the per- 
turbation expansion with respect to small A, the latter 
is surprising since at U — the slope is always positive 
for any A > 0.^'^ 

At the MIT the behavior of the self-energy changes. 
In particular, the imaginary part becomes vanishingly 
small at ti> = whereas the real part is finite. These re- 
sults imply that the mechanism for opening a correlation 
(Mott) gap at this MIT transition is different from that in 
the pure Hubbard model with particle-hole symmetryii^ 
Namely, consider the spectral function at the Fermi level, 
which is expressed in terms of the self-energy: 



A{0) 



deNoie) 



ImS(O) 



[e - RcS(0)]2 + [ImS(0)]2 



(10) 



Since we obtained numerically that at the Fermi level 
the imaginary part of the self-energy vanishes and the 
real part is larger then the band-width W, we find 
A{0) = iVo[Rci;(0)] = 0. This result is in contrast to 
the Mott-Hubbard MIT in the pure Hubbard model at 
half-filling with particle-hole symmetry. In this last case 
the opening of a correlation gap at T = is due to the for- 
mation of a delta-like singularity in the imaginary part 
of the self-energy at the Fermi level. The real part of 
S(ci;) has a 1/w divergence at this point, consistent with 
the Kramers-Kronig relations. In the disordered case we 
do not see such a behavior of RcS(lj), which implies that 
lifting of the particle-hole symmetry due to finite disor- 
der and noninteger filling affects the way how the gap is 
opened in the Mott insulator. 



V. ASYMPTOTIC LIMIT A ^ oo 

Our understanding of the Hubbard model with binary 
alloy disorder is based on the fact that the lower and up- 
per alloy subbands are split at large A. We now show 
that for A — > oo the upper (or lower) alloy subband can 
be neglected and the problem can be mapped onto a low- 
energy subspace of the full Hilbert space. Effectively, the 
correlated and binary-alloy disordered electronic system 
is represented by a correlated pure system with renor- 
malized parameters. 



A. Mapping of the Hilbert space 

We consider the case A ^ max{U,W). Then the 
Hilbert space can be divided into two subspaces A and B 
consisting of ions with energies ca and cb, respectively. 
We denote the projection operator onto the A-subspace 
by P and the projection operator onto the B-subspace 
by Q = 1 — The Schrodinger equation can be 
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FIG. 5: (a) Spectral function, (b) imaginary part of the self- 
energy, and (c) real part of the self-energy for the Hubbard 
model at V = 0.5, U = 3 and different disorder A. As A 
increases, the quasiparticle peak appears and then vanishes, 
signaling a transition from a metallic to an insulating phase. 
In the insulating phase ( A = 2) a Mott gap is opened at w = 
with ImE(O) = whereas ReE(O) remains finite. 



decomposed as 

{E ~ PHP)P\^/) = PHQQ\^) 
{E~QHQ)Q\^) ^QHPP\^), (11) 

where j^*) is a many-body eigenstate of the Hamiltonian 
with eigenvalue E. Solving this set of equations, the 
effective Hamiltonian of the lower alloy subband is for- 
mally given by 



H,s - PHP - PHQ 



QHQ 



QHP. 



(12) 



The mapping implies that the number of lattice sites cor- 
responding to the A-subspace is equal to the number of 
tA ions, i.e., Na — xNl- Therefore, if we restrict our- 
selves to this low-energy many-body subspace -P|^), and 
consider only the effective Hamiltonian iJoff, the filling 
factor can be replaced by ly* = v/x. 

The interaction U is not changed by the projection of 
the Hilbert space because it is a local quantity. However, 
the band-width of the lower alloy subband is renormal- 
ized because the number of nearest neighbors with onsite 
energy ca is reduced. In order to estimate how the band- 
width renormalizes we consider the Bethe lattice with a 
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FIG. 6: (a) Spectral function, (b) imaginary part of the self- 
energy, and (c) real part of the self-energy for the Hubbard 
model at A = 2 and U — 3 and 2. 

finite coordination number z and calculate the second 
moment /i^-* = (AjiJQ |A) for the non-interacting Hamil- 
tonian Hq ~ ^ X^i a with nearest neighbor hop- 
ping. For a given lattice site A we can have za neighbors 
with energies e^, where < za < z. It is easy to show 
that the second moment at site A is ^^^\za) — ZAt^- 
The probability distribution of za is given by 

V{zA)^(y'^^x^-{\-xY-^-. (13) 
Therefore we find the average moment 

Na 

('^^'^) = ]^E^l(^^)^(^-4) = zt^x. (14) 
i— 1 

For a Bethe lattice with coordination number z tending 
to infinity we rescale t ^ t* = tj^. We then find 
that the second moment is (/i'-^-') = V'^x. It means that 
the band- width, as measured by the second moment, is 
reduced in the effective Hamiltonian i?off by the factor 
\fx. 

B. Critical interaction Uc 

The approximate value of the critical interaction for 
the occurrence of the Mott transition can be found an- 



FIG. 7: (a) Spectral function, (b) imaginary part of the self- 
energy, and (c) real part of the self-energy for the Hubbard 
model at [/ = 2 and different disorder A. Upon increasing 
A the alloy subbands are split and a quasiparticle resonance 
emerges between the Hubbard subbands. 



alytically within the linearized DMFT, where the full 
DMFT problem is mapped onto the two-site SIAM 
with self-consistency conditions^ For the pure Hubbard 
model it was s hown that the critical interaction has the 
value Uc — Q\/JiJ^. From the results in the last subsec- 
tion we find that the critical interaction for the MIT in 
the strong disorder limit A — > oo is given by Uc = 6t*^/x. 
For parameters used in our calculations t* ~ 0.5 and 
X = 0.5 we obtain Uc = 3/V2. This value is shown by 
the dotted line in Fig. Q). The agreement between our 
numerical calculation and this estimate of Uc is surpris- 
ingly good. 



VI. CONCLUSIONS 

In this paper we studied the Mott-Hubbard transition 
in a correlated electronic system with binary alloy dis- 
order. By numerically solving the DMFT equations at 
T = we showed that for filling factor v — x a. MIT 
occurs when the disorder A or the interaction U are in- 
creased. Regarding the excitation spectrum for the elec- 
trons we introduced the notion of the "alloy Mott insu- 
lator" for J7 < A, and of the "alloy charge transfer insu- 
lator" for [/ > A. This classification is analogous to the 
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Zaanen-Sawatzky-AUen scheme for two band-Hubbard 
systems]^ In our case however, the role of the oxygen 
band is played by the upper alloy subband. We also found 
hysteresis upon approaching the metal-insulator bound- 
ary, depending on the initial conditions imposed in the 
iterative solution of the DMFT equations. It shows that 
hysteresis is a generic feature of the MIT in pure and in 
binary alloy disordered systems within the DMFT sce- 
nario. We also found that the opening of a Mott gap 
is associated with the disappearance of the imaginary 
part of the self-energy at the Fermi level. Finally we 
discussed the analytical theory, valid in the alloy band- 
splitting limit, and showed that the Hubbard model with 
binary alloy disorder can be mapped onto an effective 
Hubbard Hamiltonian with renormalized bandwidth and 
filling factor. The estimated critical interaction Uc in this 
asymptotic theory agrees very well with the numerical re- 
sults. 

Our study of the Mott-Hubbard MIT was limited to 
the Hubbard model with a non-degenerate orbital. A 
similar phase transition should be expected in the Hub- 
bard model with orbital degeneracy An important ques- 
tion is, however, whether one can find a physical system 
where such a transition is realized. For binary alloys this 
might be a very demanding task because the predicted 
MIT requires fine tuning of the filling factor with concen- 
tration of the alloy elements as well as special values for 
the interaction and the disorder splitting. At present we 
do not know which alloy system would be the best can- 
didate for realizing the predicted Mott-Hubbard MIT. 

The most promising candidates for experimental re- 
alization may come from systems with fermions moving 
on artificial lattices. Creating a lattice with a binary al- 
loy disorder seem to be possible either with a matrix of 
quantum dots with two different sizeS)^ or with optical 
lattices, where counterpropagating laser beams can be 
used to trap fermionic atoms In the latter case, with 
proper selection of laser light and physical boundaries one 
can obtain at least quasi-disordered systems^ where a 
binary alloy-like distribution is possible. From the point 
of view of our theory, where Anderson localization is not 
included, quasi-randomness of the optical lattice is not 
a major limitation as it yields the effective alloy band 
splitting which is crucial for our calculation. Since in the 
present Mott insulator the long-range ordering is sup- 
pressed due to the disorder, such a system might be very 
useful to study the ground state of, and excited states in, 
paramagnetic Mott insulators. 
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APPENDIX: CALCULATION OF THE 
SELF-ENERGY 

The self-energy can, in principle, be calculated directly 
from a Dyson equation: E(w) = G~"'^(w)— w-l-ry(w). How- 
ever, Bulla et aZf^ have proposed a different approach, 
namely to use a two-particle correlation function to de- 
termine the self-energy. In this way certain systematic 
errors, leading to inaccurate spectral weights, are can- 
celled. In this Appendix we show how to generalize this 
method to find the self-energy in the correlated system 
with binary alloy disorder. With a new interpretation of 
the CPA equation one can use two-particle correlation 
functions to determine the single-particle self-energy. 
Our method can be applied to an arbitrary DOS and, 
as we checked independently, leads to better convergence 
of the DMFT equations due to apparent error cancella- 
tions. 

Within the DMFT, the disordered Hubbard model is 
mapped onto a single-impurity Anderson model which 
contains the local ionic energy as a parameter. For 
each ei the Anderson model is solved independently. It 
yields the local e^-dependent Green function G{uj, e^) = 
{{aia\al^)). In addition, for each we can find a two- 
particle Green function 



(A.I) 



On the other hand, the equation of motion for G{Lu,ei) 
obtained from the ei-dependent SIAM is 

{iu - e,)G{Lu, - UF{u;, e,) - 7?(c^)G(c^, e,) = 1, (A.2) 

which can be rewritten as 



[to - 77(w)]G(c^, €,) = 1 + V^iu;)Giu;, e,), 



(A.3) 



where we define a complex and frequency dependent (dy- 
namical) scattering potential 



U 



G(^,e.)' 



(A.4) 



From a formal point of view equation ljA.31) looks like 
a single-particle equation for a Green function G{Lu,ei) 
with random potential Vi^Lu). The average Green func- 
tion G{lo) is found within the CPA by demanding that 
the average of the transfer matrix, given by 
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T,(lo) 



S(^) 



1 



+ '-^ G{uj,ei) 



S(w) G(w 



(A.5) 



with self-energy S(u;), vanishes, i.e., J d€iV{ei)Ti{uj) 
0. 
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For binary alloy disorder the equation for the self- 
energy can be solved analytically, leading to 



where 



Vl{Lu)+V2iuj) 



1 



(A.6) 



Vliu)+V2{LU) 



G{co) 



Vi{lu)V2{uj) ~ X 



(1 



r 



Causality of the Green function requires that the sign in 
Eq. (|A.6|1 has to be properly chosen: (i) at — oo 

(+oo) the physical solution of ljA.6|) has a (— ) [(+)] sign; 



(ii) the change of the sign happens an odd number of 
times at frequencies ojo for which ImS(wo) = 0. 
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